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By applying a variational principle on a magnetic system within the framework of extended ir¬ 
reversible thermodynamics, we find that the presence of a temperature gradient in a ferromagnet 
leads to a generalisation of the Landau-Lifshitz equation with an additional magnetic induction 
held proportional to the temperature gradient. This held modulates the damping of the magnetic 
excitation. It can increase or decrease the damping, depending on the orientation of the magnetisa¬ 
tion wave-vector with respect to the temperature gradient. This variational approach conhrms the 
existence of the Magnetic Seebeck ehect which was derived from thermodynamics and provides a 
quantitative estimate of the strength of this effect. 


I. INTRODUCTION 

The effect of a thermal spin torque on the magnetisa¬ 
tion dynamics has attracted a lot of attention recentljDll^. 
In a conductor, the spin dependence of the transport 
properties implies that a heat current induces a spin c ur- 
rent, and consequently, a torque on the magnetisation.l^ 
In an insulator, this transport model does not apply. 

In this publication, we show that extended irreversible 
thermodynamics leads to a variational principle for the 
magnetisation which predicts the existence of an addi¬ 
tional magnetic induction field proportional to the tem¬ 
perature gradient in the Landau-Lifshitz equation. In our 
previous woriP, we called this effect the “Magnetic See¬ 
beck effect” since the Seebeck effect refers to the presence 
of an electric field induced by a temperature gradient. 
This effect should not be confused with the transport 
phenomenon known as the spin Seebeck effeciPt^i^. 

Classical irreversible thermodynamics (CIT^pHU 

re¬ 
quires the system of interest to be at local equilibrium. 
Transport phenomena are then described by phenomeno¬ 
logical relationships between current densities and gen¬ 
eralised forces so as to fulfill the second law of thermo¬ 
dynamics. When a system does not satisfy the con¬ 
dition of local equilibrium, it can be described within 
the framework of extended irreversible thermodynamics 
(EIT) where the current densities are considered as ad¬ 
ditional state variableJi^. In this article, we show that 
this approach provides an expression for the Magnetic 
Seebeck effect in terms of the thermal properties of the 
magnetisation. 


II. VARIATION OF THE INTERNAL ENERGY 

In the absence of a magnetic excitation field, the 
magnetisation M is collinear to the magnetic induction 
field obtained by performing the variation of the inter¬ 
nal energy with respect to the magnetisation Su/SM, as 
pointed out by Gurevich.^^ In the presence of a magnetic 
excitation field b, the Landau-Lifshitz equation describes 
the precession of the magnetisation M about this mag¬ 
netic induction field. Since the magnetisation is locally 


out of equilibrium, we use the framework of extended irre¬ 
versible thermodynamics. According to this framework, 
the internal energy density u {M, V x M) is a function 
of the magnetisation M and of the magnetisation current 
j M — ^ X M that are in turn functions of the position 
r. According to the variational equation (A9) established 


explicitly in the Appendix, the variational derivative of 
the internal energy in the bulk of the system reads. 


6u 

SM 


du 

Mi 


-I- V X 


du 


a(v X M) 


( 1 ) 


The variational principle used by Gurevich et and 
Bose et am assumes that the internal energy density is 
a function of the magnetisation M and the gradient of 
the magnetisation V M. The choice of the curl of the 
magnetisation V x M as a degree of the freedom of the 
internal energy density is motivated by the framework of 
extended irreversible thermodynamics. Furthermore, for 
a system where the magnetisation is driven out of local 
equilibrium, the internal energy density is expressed in 
the bulk as,^iSl 


u (M, V X M) = - (V X M) • A (M) 


( 2 ) 


where A (M) is the potential vector. This implies that 
the second term on the RHS of the variational equa¬ 
tion Q corresponds to a magnetic induction field, as 
it should. 

We consider a slab subjected to a temperature gradi¬ 
ent V T and a uniform and constant magnetic induction 
field Bq that are applied in plane along the i-axis as 
illustrated in Fig. The magnetic excitation field b is 
applied orthogonally to the i-axis. The constant applied 
magnetic induction field Bq and the magnetic excitation 
field b are oriented as follows. 


Bq ■ Z — Bq 
bBo = 0 


( 3 ) 


The magnetisation M is the sum of the saturation 
magnetisation Ms and of the magnetic response field 

m, i.e. 


M = Ms + rn 


( 4 ) 
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where the norms satisfy the relations 

\\Ms\\ » ||m|| 

IITWII = Mg + = cste 

and their orientations are given by, 

M s ■ z = Mg 
m ■ Mg = 0 


( 5 ) 


( 6 ) 


The term on the LHS and the first term on the RHS 
of equation Q are recast as, 


5u 

Jm 

du 

5M 


— ^0 + 

= Bq + b 


( 7 ) 


where the first-order magnetic induction field Si is ori¬ 
ented as follows, 


Bi ■ Bq = 0 


( 8 ) 


The temperature gradient V T is imposed along the z- 
axis as illustrated in Fig. [T] Along that axis, the spatial 
symmetry of the system is broken, i.e. 


VT - 9 

V = z — 
dz 

The conditions ([^ and imply that, 
V X M = V X m 


( 9 ) 


( 10 ) 


The relations Q and (10) imply that the variational 
equation 0 is recast as. 


Si = 6 -I- V X 


du 

d (V X m) 


( 11 ) 


III. MAGNETIC SEEBECK EFFECT 


The second term on the RHS of the relation (11) is 
recast as, 


Vx 


du 


9 (V X m) 


d{V X m) 


= I I 11 


-1 


In the linear response, i.e. to first-order in the mag¬ 
netic excitation field m, the differential operator d/du 
commutes with the differential operator V which implies 
that, 


d (V X m) dm 

— - - = V X - 

du du 


(13) 


Substituting the relation (131 into the expression (12) 
yields. 


V X 


du 


Vx(V-^x|^) (14) 

9 (V X m) J \ dm J 


In the linear response, the relations 0-0 imply that 
du 


dm 


= b 


Taking into account the conditions 0, (15) and 

V • = 1 


the expression (|14|) becomes 

V X 


du 


d (V X m) 


= -i-v-i ( —b 

dz 


d 


which implies that the relation (11) is recast as, 

d 


B, 


b- z-V-^{^b 
oz 


The linear magnetic constitutive equation reads, 

b = 


which implies that the relation (|18|) becomes, 
Bi=b- 




The derivative of the second condition 0 yields. 


dMg dm 

Mg + m — = 0 
dz dz 


(15) 


(16) 


(17) 


(18) 


(19) 


( 20 ) 


( 21 ) 


which implies that for a uniform precession in the plane 
orthogonal to the £-axis. 




dz 


\ dz J 


m 


( 22 ) 


The saturation magnetisation Mg = Mg (T) is a function 
of the temperature, which implies that. 


dMg dMg 


dz 


dT 


(z-VT) 


(23) 


The magnetic susceptibility x is a function of the satu¬ 
ration magnetisation Mg, i.e. 


dx ^ dx dMg 
dz ~ dMg dT 


(z-VT) 


(24) 


In the linear response, the relations (22)-(24) imply 
that the expression (20) is recast as, 


= 6 


^ / 1 ^ / dMg^\ , 

xKxdMg^ j{ dT > 


(25) 

For a ferromagnet, the magnetic susceptibilit y y is pro¬ 
portional to the saturation magnetisation MJ^. Thus, 
the first condition 0 imposes on the first term in brack¬ 
ets (25) the condition, 


1 

X dMg Mg 


(26) 
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Taking into account the condition (26), the relation (25) 
reduces to, 


and ( [M| ) yields the following expression for the Magnetic 
Seebeck parameteil^. 


Bi = b 


Mo Ms f dMs \ 
Xm? \ dT ) 


(VT) • 


(27) 


MoM| dM^\ 0 

nks \MsJ \xMs dT ) 


(31) 


The first-order magnetic induction field Bi that is or¬ 
thogonal to the zeroth-order magnetisation saturation 
Ms exerts a thermal magnetic torque r = Ms x B^ 
on the magnetisation M illustrated in Fig. 



FIG. 1: Precession cone of the magnetisation M where the 
saturation magnetisation Ms, the constant magnetic induc¬ 
tion field Bo and the temperature gradient V T are oriented 
in plane along the z-axis. The magnetic response field m and 
the magnetic torque r are orthogonal to the z-axis. 


To compare the result of this analysis with our previ¬ 
ous work on the Magnetic Seebeck effeclP, we recast the 
relation (27) as. 


Bi — b — fj,o {kr ■ V 


(28) 


where the thermal wave-vector is given by, 

feT = -^(^)VT (29) 


= _ils,(^)vr 

Xm^ \ dl J 


where dMs/dT < 0. The expression (28) for the mag¬ 


netic induction field Bi has the same structure as equa¬ 
tion (7) of our previous worlPI. However, in the second 
term on the RHS accounting for the Magnetic Seebeck 
effect, the expression of the thermal wave-vector dif¬ 
fers. According to our previous work in the classical ir¬ 
reversible thermodynamical frameworlP, the expression 
for the thermal wave-vector is. 


A n fcs ^ ^ 

KT = -7T9 V 1 

Mo 


(30) 


where n is the Bohr magneton number density, ks is 
Boltzmann’s constant and A is the dimensionless Mag¬ 
netic Seebeck parameter. Identifying the relations (29) 


since dMs/dT < 0. The thermal vector kx generates a 
magnetic induction field (28) that is proportional to the 
temperature gradient V T. This field leads to a general¬ 
isation of the Landau-Lifshitz equation. The respective 
orientation between the wave-vector k of the magneti¬ 
sation waves and the temperature gradient V T leads 
to an increase or attenuation of the magnetic damp¬ 
ing. As predicted theoretically and observed experimen¬ 
tally for magnetostatic backward volume modes in a YIG 
slalP, the magnetic damping is attenuated for magnetisa¬ 
tion waves propagating along the temperature gradient 
while it is increased for magnetisation waves propagat¬ 
ing against the temperature gradient. This modulation 
of the magnetic damping is a consequence of the Mag¬ 
netic Seebeck effect. 


IV. CONCLUSION 

In this article, the variational approach for the descrip¬ 
tion of magnetisation dynamics is performed in the ex¬ 
tended irreversible thermodynamical framework. It pre¬ 
dicts the existence of a Magnetic Seebeck effect for the 
propagation of magnetisation waves along a temperature 
gradient in a magnetic slab. In the classical irreversible 
thermodynamical framework, the coupling between the 
magnetisation dynamics and the thermal gradient is ex¬ 
pressed by a phenomenological relation imposed in or¬ 
der to satisfy the second principle of thermodynamics. 
By contrast, in the extended irreversible thermodynam¬ 
ical framework, the coupling between the magnetisation 
dynamics and the thermal gradient is intrinsic to the de¬ 
scription of the system itself, i.e. it derives from the ther¬ 
mal properties of the system. The comparison between 
the expressions obtained for the thermal wave-vector kx 
in both frameworks yields an explicit expression for the 
dimensionless parameter A that defines the strength of 
the Magnetic Seebeck effect. 


Appendix A: Variational principle 

In a stationary state, the kinetic energy associated 
with the precession is constant. Thus, it can be 
ignored while performing the action variation. The 
action S {M, V x M) and the Lagrangian density 
C {M, V X M) are functions of the magnetisation M 
and of the magnetisation current V x M that are in 
turn functions of the position r. These quantities are 
related by the integral expression, 

S {M,V X M) = [ dtd^rC{M,V X M) (Al) 
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In a stationary state, the internal energy density 
u {M, V X M) is equal to the opposite of the Lagrangian 
density C {M, V x M) up to a constant, i.e. 


C (M, V X M) = -u (M, V X M) 


(A2) 


where the internal energy density plays the role of the 
potential. Thus, the action (Al) is recast as. 


S{M,V X M) = - J dtd^ru{M,V X M) (A3) 

The variation of the action S {M, V x M) is expressed 
formally as. 


SS = — dtd^rSu 


and yields. 


6S=- dtd-^ 


du 


■ 6M 


du 


S{V X M) 
(A5) 


dM ' d(y X M) 

Using the vectorial identity, 

A ■ S {V X M) = {V X A) ■ SM - V ■ {A X SM) (A6) 


the relation (A5) is recast as 
du 


6S = - dtd-^r 


dM 


+ / dtd^rV 


+ V X 
du 


du 


d{V X M) 


9(V X M) 


X 6M 


SM 

(A7) 


Using the divergence theorem, the second integral on the 
RHS of the relation (A71 is recast as a surface integral. 


Taking the limit where the integration surface tends to 
infinity and assuming that the magnetisation is uniform 
at infinity, this integral vanishes, which implies that the 
relation (A7) reduces to. 


6S = - 


dt d^r 


du 

dm 


+ V X 


du 


a(v X M) 


(A4) Identifying the integrands in relations (A4) and (A8) 


• SM 
(AS) 


and taking the variational derivative of the functional 
u {M, V X M) with respect to Af, also known as func¬ 
tional derivativ^^ yields, 


Su 

SM 


du 

m 


-I- V X 


du 


5(V X M) 


(A9) 
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